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ABSTRACT 

Context. High-resolution studies of class 0 protostars represent the key to constraining protostar formation models. VLA16234-2417 
represents the prototype of class 0 protostars, and it has been recently identified as a triple non-coeval system. 

Aims. We aim at deriving the physical properties of the jets in VLA16234-2417 using tracers of shocked gas. 

Methods. ALMA Cycle 0 Early Science observations of CO(2-l) in the extended configuration are presented in comparison with 
previous SMA CO(3-2) and Herschel- PACS [Oi] 63 pm observations. Gas morphology and kinematics were analysed to constrain 
the physical structure and origin of the protostellar outflows. 

Results. We reveal a collimated jet component associated with the [Oi] 63 pm emission at about 8" (~960 AU) from source B. This 
newly detected jet component is inversely oriented with respect to the large-scale outflow driven by source A, and it is aligned with 
compact and fast jet emission very close to source B (about 073) rather than with the direction perpendicular to the A disk. We also 
detect a cavity-like structure at low projected velocities, which surrounds the [Oi] 63 pm emission and is possibly associated with the 
outflow driven by source A. Finally, no compact outflow emission is associated with source W. 

Conclusions. Our high-resolution ALMA observations seem to suggest there is a fast and collimated jet component associated with 
source B. This scenario would confirm that source B is younger than A, that it is in a very early stage of evolution, and that it drives a 
faster, more collimated, and more compact jet with respect to the large-scale slower outflow driven by A. However, a different scenario 
of a precessing jet driven by A cannot be firmly excluded from the present observations. 
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1. Introduction 

With lifetimes less than 10 5 yr, class 0 objects represent the ear¬ 
liest phase of star formation, whe n most of their mass is still in 
the fo rm of dense envelopes (e.g. lAndri et al. 2000; E vans et ail 
2009). Because they have not undergone significant evolution, 
they are likely to retain information on the physical and chem¬ 
ical initial conditions and on the physics of the collapse phase. 
High-resolution studies of class 0 protostars thus represent the 
key to constraining protostar formation models. 

Located in p Ophiuchus (d~ 120 nc, iLoinard et alJ 12008 ), 
VLA16234-2417 (hereafter V LA1623) is considered to be the 
prototype of class 0 protostars dAndre et al.ll 199(1119931) . It was 
originally identified as a strong continuum source at the mm 
wavelengths associated w ith a large-scale bipolar outflow tra ced 
by CO and Hi emissions ([Andre et al.il9901:iDent et al.l 1995b Yu 
& Chernin 1997) nearly on the plane of the sky (1=75°, Davis et 


multiple protostellar system (see also [Maury et al. 201 Of [201 2j). 
H owever, VLA1623N multiplicity ([Looney et al.l 20001: Dent et 
al. 1 19951: Mi&Chemin|l997) was recently confirmed bv Murillo 
& Lai (1201 3 all and lChen et al.| (|2013[) through SMA observa- 
tions. In particular. iMurillo & Lail (120 1 3al) show that VLA1623 
represents a triple non-coeval system composed of VLA1623A, 
VLA1623B, and VLA1623W (hereafter A, B, and W), with each 
source driving its own outflow. Source A was identified as a 
deeply embedded class 0 source with no emission at wavelengths 
shorter than 24 pm. It shows a rotationally supported disk in 
C ls O and C 17 0, later confirmed through ALMA observations by 
iMurillo et al. ( 2013b ). and drives the large-scale bipolar outflow. 
It has been suggested that source B is a cold and compact source 
in a very early stage of star formation between the starless core 
and class 0 stages with a separation of l'/l (~ 1 3 2 AU) from A 
see also Loo nev et al.ll2000t iMaurv et al 1l20 1 2l) . IMurillo & Lail 


al. 19991) . High-resolution observations bv [Bontemps & Andrei 
(1 1997b revealed multiple continuum sources, which were inter¬ 
preted as a radio jet driving the large-scale outflow rather than a 


[201 3al) suggest that source B drives a pole-on outflow. Finally, 
it has been argued that source W is an L-IA L e class I source at 
the projected distance of 10" from A. 
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Fig. 1 . ALMA C0(2-l) emission integrated over the blue- (-8, +1 km s *) and red-shifted (+5, +18 km s *) velocity ranges are compared with 
the PACS [Oi] emission line at 63 pm bv lNisini et alJ ( 2013b . The ALMA synthesised beam of the CO map is shown in the right panel. Green 
triangles indicate the positions of the three sources, as reported by Murillo et al. (201 3bj) ; magenta triangles label the closest CO kn ots identified by 
Andr e et aTI d!99Ch . Black lines indicate the directions of the VLA1623A disk fPA=3.T. lMurillo & Lail2013HlMurillo et al.l20 1 3bll and the relative 
outflow (PA=125°). The latter is consistent with the large-scale CO outflow direction and with the secondary [Oi] 63 pm peaks, whereas the closest 
[Oi] 63 pm emission peaks are misaligned. The cyan line indicates the direction of the collimated jet (Brl and Bbl CO knots from Fig. EK which 
is not associated with the large-scale CO outflow direction and with the direction perpendicular to the A disk, but it is rather associated with the 
closest [Oi] 63 pm emission peaks. 


In this research note, we present Atacama Large Millime¬ 
ter/submillimeter Array (ALMA) Early Science Cycle 0 obser¬ 
vations of CO(2-l) towards VLA1623. The goal is to study 
the physical properties of the jet emission in VLA1623 at high- 
angular resolution. 


2. Observations 

VLA1623 has been mapped in the CO(2-l) line with ALMA 
during the Early Science Cycle 0 phase on April 8, 2012 (PI 
Murillo). The details of the observations, calibra tion, and data 
reduction were presented in lMurillo et alJ ( 2013bh . Pointing co¬ 
ordinates were <+ 12000 = 16 h 26 ln 26'.419 and <5 12000 =—24°24 / 29799. 
Observations were done in Band 6 (230 GHz) using the extended 
configuration, and they consisted of 16 antennae with a baseline 
range of ~36—400 m for a total observing time of one hour. The 
minimum baseline of 36 m yields a maximum recoverable scale 
of 4" at the selected frequency, which is consistent with filtering 
most of the emission associated with the large-scale outflow (see 
Sect. [4]). The synthesised beam is 0774x0758 (PA=93°). The rms 
noise is 3-5 mJy beam -1 per 1 km s -1 channel. Data was anal¬ 
ysed using the GILDASQ package. 


3. Results 

Figure [T] presents the comparison between the ALMA blue- and 
red-shifted CO(2-l) emission, the large-scale CO outflow emis¬ 
sion, and the PACS [Oi] 63 pm line emission from lNisini et al.l 
d2015l) . The [Oi] 63 pm emission is extended along the outflow 
direction and shows several emission peaks corresponding to dif¬ 
ferent shock episodes. The direction perpendicular to the A disk 


1 http://www.iram.fr/IRAMFR/GILDAS 


(black solid line, PA= 125°) is consistent with the large-scale out¬ 
flow direction (identified in Fig. (T|by the BLUE1/RED1 peaks) 
and with the two associated secondary [Oi] 63 pm peaks. Con¬ 
versely, the closest [Oi] 63 pm peaks are not co-aligned with the 
A outflow direction, but they are instead aligned with source B 
(cyan line). 

The ALMA data resolve the CO emission previously de¬ 
tected in CO(3-2) with SMA (see lMurillo & Laill2013a[ and 
Fig. IA.2I where the channel maps of the SMA observations are 
presented); in particular, the southern (low-velocity) blue-shifted 
emission is almost completely filtered out by the ALMA obser¬ 
vations. Moreover, no compact outflow emission seems to be 
associated with source W. Finally, the ALMA CO(2-l) emis¬ 
sion is not associated with the large-scale outflow emission, and 
it possibly shows two outflow components: a collimated emis¬ 
sion, identified by the northern blue- and southern red-shifted 
emissions along the cyan solid line, which is inversely oriented 
with respect to the large-scale outflow and associated with the 
closest [Oi] 63 pm emission peaks; and an arc-like feature of 
enhanced red-shifted CO intensity, which surrounds the north¬ 
ern [Oi] 63 pm emission and may be associated with a cavity 
opened by the jet. 

More clues about these components are provided by the 
channel maps of the red- and blue-shifted CO(2-l) emission 
presented in Figs. [2] and IA.1I The two components identified 
in the velocity-integrated emission (Fig. Q]) appear to be well 
separated in velocity, the red-shifted cavity-like structure dom¬ 
inating the emission at low projected velocities (Arl, Ar2, and 
Ar3 knots), whereas the collimated jet associated with the clos¬ 
est [Oi] 63 pm peaks emitting in the higher-velocity channels 
(Br2, Brl, Bbl, and Bb2 knots). 

Finally, compact blue- and red-shifted outflow emission at 
about 073 from B are detected (Bbl and Brl knots). We note that, 
assuming an inclination angle with respect to the line of sight 
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Fig. 2. Channel maps of the ALMA C0(2-l) emis sion integrated over three velocity ranges: systemic velocity in the left panels (v sys ±l km s 1 , 
where the systemic velocity is 3.6 km s -1 , see lYu & Chemin|l 19971 [Narayanan & Loganll2006l) . low-velocity in the middle panels (LV; -2, +3 
km s~* and +5, +10 km s~*), and high-velocity in the right panels (HV; -8, -2 km s -1 , and +10, +18 km s~*). Offsets are given with respect to 
the pointing coordinates, i.e. aj2ooo=16 h 26 m 26)419 and (5j2ooo=-24 0 24'29'.'99. Contour levels of emission are traced at 3, 5, and 8 cr levels and 
increase in steps of 5 cr, with the exception of the HV panels where they are shown at 3 and 5 cr and increase in steps of 10 cr. Negative contours 
are displayed with dotted (upper) and dashed (lower) lines starting at 3 cr level and decreasing in steps of 3 cr. The ALMA synthesised beam of 
the CO map is shown in the left panels. Black triangles represent the three protostellar sources, which are labelled in the left panel. Other labels 
indicate the CO features discussed in the text. Black and cyan solid lines are the same as in Fig.[7| 


of 75° as for the A outflow (iDavis et al.lll999 ). the maximum 
projected velocities of 14 km s _l (with respect_to the ambient 
systemic velocity of 3.6 km s~L Yu & Cherninl 19971: Narayanan 


& Logan 12006I) . detected in association with the collimated jet 
component, correspond to de-projected velocities of 54 km s -1 . 
Even assuming an inclination angle of 45°, de-projected veloc¬ 
ities of about 20 km s _1 are derived. Thus, we refer to the col¬ 
limated jet component associated with the closest [Oi] 63 /urn 
peaks and to the compact jet emission from B as fast jet compo¬ 
nents. 


4. Discussion 

We observed compact jet emission very close to B (Bbl and Brl 
knots). In addition, farther from sources A and B we detect two 
outflow components: one faster red- and blue-shifted jet emis¬ 
sion associated with the closest [Oi] 63 fim peaks and inversely 
oriented with respect to the large-scale outflow emission (Bb2 
and Br2), and a slower northern red-shifted cavity-like structure 
surrounding the [Oi] 63 gm emission (Arl, Ar2, and Ar3). Two 
possible scenarios may explain the present observational results: 
1) precession of the A jet, and 2) the presence of a fast compact 
jet associated with source B. 

In the first case, the misalignment of the closest [Oi] 63 n m 
emission peaks with respect to the direction of the large-scale 


outflow may be explained with a scenario in which the A jet, 
with an inclination angle close to the plane of the sky, is precess- 
ing owing to the mutual interaction with the B companion. This 
causes a change in propagation direction of the outflow driven 
by A and the inversion of the blue- and red-shifted emissions. 
We note that the presence of precessing jets in this region has 
already been suggested bv ICaratti o Garatti et alJ (2 006h . but on 
much larger scales than those probed by the present ALMA data. 
In this view, the compact jet emission associated with the ex¬ 
tremely young source B (Bbl and Brl knots)_would represent a 
pole-on outflow, as already suggested bv lMurillo & Lain2013 a). 
However, the direction of the secondary, hence (in this view) of 
older [Oi] 63 gm peaks, is consistent with the direction perpen¬ 
dicular to the A disk (and with the large-scale outflow), whereas 
if the A jet was precessing we would expect the direction of the 
closest [Oi] 63 yum peaks to be consistent with the direction of 
the A jet. 

A second interpretation could be that the fast collimated jet 
emission associated with the closest [Oi] 63 /rm peaks (Bb2 and 
Br2) and the high-velocity knots close to source B (Bbl and Brl) 
represent a newly detected jet driven by B. This is supported by 
the position-velocity diagram along the direction of this fast col¬ 
limated jet (identified by the Brl and Bbl knots), which is pre¬ 
sented in Fig. [3] It shows how the pairs of corresponding red- and 
blue-shifted knots lying along this possibly new jet seem rather 


Article number, page 3 of[8] 

































A&A proofs: manuscript no. vla!623alma 



-10 -50 5 10 

Offset (arcsec) 


Fig. 3. Position-velocity cut of the CO(2-l) emission in false colours and black contours performed along the direction of the collimated jet 
identified by Brl and Bbl CO knots in Fig.[2](along the cyan solid line in Figs. EH21 PA=116°). Labels are the same as in Fig. [2] Dashed vertical 
white lines mark the positions_of VLA1623A and VLA1623B, while the dashed horizontal white line indicates the ambient systemic velocity of 
3.6 km s ~‘, d erived bv|Nara yanan & Logan ( 2006 1 from N 2 H + observations taken with the Five College Radio Astronomy Observatory (FCRAO) 
(see also lYu & Cherninll 19971 ). Dotted black lines connect the corresponding blue- and red-shifted knots associated with the possibly new fast and 
collimated jet emission detected with ALMA and labelled in Fig. [2] Both pairs of corresponding knots seem to be associated with B rather than 
with A. 


to be associated with source B than with source A. This would 
mean that the jet_associated with B is not pole-on as previously 
argued by iMurillo & Lail d2013al) . and it implies a very short 
dynamical timescale for the close compact emission knots (Brl 
and Bbl). In particular, assuming an inclination angle of 75° as 
for the A outflow dDavis et alil!999f) , the dynamical timescale 
of Brl and Bbl, located at the projected distance from B of 073 
and moving at the radial velocity of ~ 12 km s 1 (Figs.[2]-[3]), is 
~4 yr. Even assuming an inclination angle of 45° for the B jet, 
the dynamical timescale would be ~14 yr. 

On the other hand, the slower cavity-like structure would be 
associated with the large-scale outflow driven by A, which is 
possibly extended and poorly collimated and therefore almost 
completely filtered out by the ALMA observations. We note that 
the lack of CO emission in our map near the systemic velocity 
is probably due to interferometric spatial filtering of extended 
emission from the envelope and the large-scale outflow, there¬ 
fore our CO map can only trace emission from compact red- 
and blue-shifted material. According to this scenario, source B 
is driving a much more compact and collimated outflow with re¬ 
spect to A. In fact, if we assume for the B jet the same intrinsic 
length as observed for the A outflow, i.e. at least 60-100 arcsec, 
the projected length of the B jet of ~8" would imply an incli¬ 
nation for the B jet of ~5°-8°, which can be excluded given the 
morphology of the jet emission and the well-separated jet lobes. 
In conclusion, this interpretation is consistent with source B be¬ 
ing younger than A, in_a_very early stage of evolution as already 
proposed bv IMurillo & Lail ( 2013 a). and driving a faster, more 
collimated, and more compact jet with respect to the large-scale 
slower outflow driven by A. 

Finally, the non-detection with ALMA of any outflow emis¬ 
sion towards source W could confirm that source W is older than 
the other two sources, A and B (iMurillo & Lail[2013al) . and thus 
has already swept out all surrounding material. Alternatively, the 


W outflow emission, if not compact, could be filtered out by the 
interferometer. 

ALMA follow-up observations of optically thin (SO) and se¬ 
lective shock tracers (SiO) would be crucial for clarifying the 
proposed scenario. In particular, observations in both the ex¬ 
tended and compact configurations are needed in order to avoid 
filtering of the large-scale emission and to eventually recover the 
emission from the large-scale outflow. This will allow the two 
possible jets to finally be disentangled and will put strong con¬ 
straints on current theoretical models of multiple system forma¬ 
tion. Moreover, wide-band spectra of the A and B sources would 
be needed to understand their nature and chemical compositions 
(from the possible detection of complex organic molecule emis¬ 
sion). 

5. Conclusions 

Compact jet emission is detected very close to source B. In addi¬ 
tion, two outflow components are detected further from sources 
A and B: one faster and more collimated jet emission associ¬ 
ated with the closest [Oi] 63 pm peaks and inversely oriented 
with respect to the large-scale outflow emission; and a slower 
cavity-like structure surrounding the [Oi] 63 pm emission. The 
observations seem to confirm that source B is younger than A, 
is in a very early stage of evolution, and is driving a faster, more 
collimated, and more compact jet, with respect to the large-scale 
slower outflow driven by A. However, a different scenario of a 
precessing jet driven by A cannot be firmly excluded from the 
present observations. 
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Appendix A: CO channel maps 

We show in Fig. IA.1I the channel maps of the blue- and red- 
shifted ALMA CO(2-l) (continuum subtracted) emission to¬ 
wards VLA1623, which are analysed and discussed in the main 
text (see Sect. [3}. The images clearly show two gas compo¬ 
nents that are well separated in velocity: a red-shifted cavity¬ 
like structure dominating the emission at low projected veloci¬ 
ties (Arl, Ar2, and Ar3 knots); and a collimated emission in the 
higher-velocity channels (Br2, Brl, Bbl, and Bb2 knots along 
the cyan line) and associated with the closest [Oi] 63 pm peaks 
(see Fig. Q}. In particular, compact high-velocity outflow emis¬ 
sion is detected close to source B (Bbl and Brl knots). The chan¬ 
nel maps of the blue- and red-shifted SMA CO(3-2) emission 
are also presented in Fig. IA.2l for comparison. We note that the 
low-velocity blue-shifted emission between 0 and +3 km s 1 is 
resolved by the ALMA observations. 
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Fig. A.l. Above: the channel maps of the ALMA CO(2-l) emission. The velocity range of each map is indicated in the top left corner. The ALMA 
synthesised beam of the CO map is shown in the central panel. Contour levels of emission are traced at 3, 5, and 10 <x levels and increase in steps 
of 10 a. Negative contours are displayed in black starting at 3 cr level and decreasing in steps of 3 <x. Black crosses represent the three protostars, 
which are labelled in the central panel. Other labels indicate the CO features discussed in the text (see also Fig. EL Black and cyan solid lines are 
the same as in Fig. |T] The ALMA observations resolve the low-velocity blue-shifted emission (0, +3 km s -1 ) detected with SMA (see Fig. IA.2l l. 
which is possibly associated with extended emission and, therefore, filtered out by our high-resolution ALMA interferometric observations. 
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Fig. A.2. Channel maps of the SMA CO(3-2) emission. The velocity range of each map is indicated in the top left corner. The SMA synthesised 
beam of 2'.'32xl'/41 (PA=55°) is shown in the central panel. Contour levels of emission are traced starting at 3 cr levels and increasing in steps 
of 3 cr. Negative contours are displayed in black starting at 3 cr level and decreasing in steps of 3 cr. Black crosses represent the three protostars, 
which are labelled in the central panel. Black and cyan solid lines are the same as in Fig.Q] 









